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Abstract 

We explore theoretically the electroluminescence of single molecules. We adopt a local-electrode 
framework that is appropriate for scanning tunneling microscopy (STM) experiments where elec- 
troluminescence originates from individual molecules of moderate size on complex substrates: Cou- 
plings between the STM tip and molecule and between the molecule and multiple substrate sites are 
treated on the same footing, as local electrodes contacting the molecule. Electron flow is modelled 
with the Lippmann-Schwinger Green function scattering technique. The evolution of the electronic 
energy levels of the molecule under bias is modelled assuming the total charge of the molecule 
to be invariant, consistent with Coulomb blockade considerations, but that the electronic occupa- 
tions of the molecular HOMO and LUMO levels vary with changing bias. The photon emission 
rate is calculated using Fermi's golden rule. We apply this theoretical approach to the STM/Zn- 
etioporphyrin/Al2O3/NiAl(110) system, and simulate various configurations of coupling strength 
between the molecule and substrate. We compare our results to the experimental observations of 
Qiu, Nazin and Ho [Science 299, 542 (2003)] for this system and find that our model provides a 
comprehensive explanation of a multitude of previously unexplained observations. These include 
the different types of current-voltage characteristics (CVC's) that are observed experimentally, the 
observed association of electroluminescence with some CVC's and not others and key properties 
of the observed photon spectra. Theoretical predictions are presented for further single-molecule 
electroluminescence experiments. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 



In the past 15 years, molecular electronics has become a field of intense interest for 
fundamental research, with potential applications in the creation of nanoscopic devices [l], 



At the same time, great progress has been made in the creation of nanoscale photonic 
devices such as those based on photonic band gap materials [4]. 

The scanning tunneling microscope (STM) is proving immensely useful in bridging the gap 



between these two fields. In STM experiments on clean surfaces, light emission has 



been observed: Systems with an STM tip over a metallic or semiconducting 



3, 



xequently 
surface 



are known to give off light due to the decay of plasmons. Enhanced photon emission has been 



|9l,[iol,[iil,Q,Q,Q. 



observed when molecules are placed inside the tip-substrate junction^ 
However, it was unclear until recently whether the stronger emission was limited to an 
amplification of the plasmon-based emission seen on metallic surfaces jisl. 16], or whether 
there could in some cases be a different, inherently molecular emission mechanism at work. 

Recently, it has been definitively demonstrated through STM experiments that electric 
current fiow through a molecule may indeed cause the molecule to luminesce [itI. [isl due to 
molecular orbital electronic transitions. This phenomenon, bridging the areas of photonics 
and molecular electronics, is a promising step towards an emerging field of single-molecule 
optoelectronics. 

Much insight into the electronic properties of these STM/molecule/substrate systems has 
been obtained by directly studying electric current, for example through the comparison of 
experimental and theoretical STM topographs and current-voltage (I-V) curves. In recent 
years, with photonic properties of STM/molecule/substrate systems also being studied, a 
photon detector and spectrometer have been added to the standard STM apparatus. Simul- 
taneous photon emission and electric current measurements have the potential to greatly 
enhance our understanding of these systems. A theoretical understanding of single molecule 



electroluminescence, however, is still in the earliest stages pj|, 120|, l2ll] and contact between 
theory and any specific experiment has not been made. The purpose of this article is to begin 
to bridge this divide between theory and single-molecule electroluminescence experiments. 

The basic idea of molecular electroluminescence as observed in STM experiments is as 
follows: By positioning an STM tip above a single molecule on a substrate, and applying 
a bias voltage between the tip and substrate, electron transmission through the molecule 
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may occur, mediated by the molecule's electronic orbitals, and the molecule may be found 
to luminesce. In a simplified picture, when a bias voltage is applied, the molecule moves 
out of equilibrium, with a fiux of electrons passing through it. If two molecular orbitals 
are located in the energy window between the electrochemical potentials of the STM tip 
and substrate, they will both be partially occupied and, if optical transitions between them 
are not forbidden, transitions from the higher-energy orbital to the lower-energy orbital will 



occur, resulting in photon emission M 

n I 

It has been predicted [19| and confirmed experimentally [18| that the relative coupling 
strengths of the molecule to the electron source and drain greatly affect molecular electro- 
luminescence. If the coupling strengths are highly asymmetric, photon emission is severely 
quenched. [3] Thus, in STM/molecule/substrate experiments where the STM- molecule cou- 
pling is normally weak, a thin insulating 'spacer' layer between the molecule and metallic 
substrate can enhance photon emission by reducing the strength of the molecule-substrate 
coupling and making it comparable with the molecule-STM coupling. Conversely, this spacer 
layer has also been shown to strongly suppress plasmon-mediated photon emission, and thus 
facilitate resolving molecular electroluminescence from the background plasmon-mediated 
photon emission that may be present even in the absence of a molecule between the STM 
tip and substrate. For instance, Qiu, Nazin and Ho have studied electroluminescence of 
Zn(II)-etioporphyrin I using a sA-thick aluminum oxide insulating layer below the molecule 
(STM/Zn-etioporphyrin/Al2O3/NiAl(110))[l3]. In these experiments, the STM image, the 
measured I-V curve and the observation of molecular electroluminescence all depend on the 
precise location of the molecule on the Al2O3/NiAl(110) substrate. 

In order to theoretically model systems such as these, where there is a thin insulat- 



ing 'spacer' layer 



27( 1 that has a complex atomic structure and a local geometry under the 



molecule that is not measured experimentally but transmits electrons nonuniformly, a local 
electrode approach has proved useful 22{]. By considering the tip-molecule and molecule- 
substrate couplings on equal footings, as /oca/ electrodes coupled to the molecule, it has been 
shown that the experimentally observed location-dependent STM images of the molecules 



can be explained in terms of different locations of dominant molecule-substrate coupling [22|. 
For the STM/Zn-etioporphyrin/Al2O3/NiAl(110) system, there is evidence that the out-of- 
plane ethyl groups of the molecule may be the locations of dominant molecule-substrate 
coupling, and that the strength of the coupling between each ethyl group of the molecule 
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and the substrate depends on the location of these groups on the substrate [22]. Thus it 
differs from molecule to molecule adsorbed on the substrate. However, to date there has 
been no theory of electroluminescence from this system. 

In this article, we extend the above local-electrode theoretical framework to the study 
of electroluminescence and I-V characteristics observed in the experimental system of 
STM/Zn-etioporphyrin/Al2O3/NiAl(110).fl^ We consider one local STM tip probe above 
the molecule, and four local substrate contacts positioned below the four ethyl groups of 
the molecule. By varying the coupling strengths between the molecule and each of the 
electrodes, differing configurations can be simulated. In this model, each electrode is repre- 
sented using a one-dimensional tight-binding model, and electron flow is modelled using the 
Lippmann-Schwinger Green-function scattering technique. Fermi's Golden Rule is used to 
calculate photon emission spectra. 

At the present time, there is no satisfactory first principles theory of the electronic struc- 
ture of molecules that are weakly coupled to the electrodes under applied bias [23|, the 
situation under consideration here: The ab initio approach to electrical conduction based on 
standard time-independent density functional theory breaks down for such systems, yielding 
unphysical behavior for the molecular energy levels and the transmission resonances asso- 
ciated with them as the applied bias is varied, and therefore incorrect calculated current- 
voltage characteristics for the molecule I23]. Thus, we ado pt a different theoretical approach: 



251] to calculate the molecular or- 



We use semi-empirical extended Hiickel parameters [2^ 
bitals and their energies at zero applied bias. The dependence of the molecular energy levels 
on the applied bias is then calculated by a self-consistent procedure based on the assumption 
that the net charge on the molecule does not change significantly as the bias applied between 
the STM tip and substrate is varied in the range of bias voltages being considered. This 
assumption is known to be appropriate for molecules weakly coupled to the electrodes, for 
example, in the Coulomb blockade regime that is not captured correctly by density func- 
tional theory. For the present system this methodology is remarkably successful, and we 
are able to attribute the prominent features of the experimental data (all of the peaks in 
the differential conductance vs. applied bias, the bias voltages at the onset of electrolumi- 
nescence and the energies of peaks in the observed photon spectra) to the movement of the 
molecular LUMO and HOMO energy levels relative to the electrochemical potentials of the 
source and drain electrodes that follows directly from the requirement that the charge on 
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the molecule is approximately independent of the applied bias. 

The experimental conductance and electroluminescence data for this system is multi- 
faceted, depending qualitatively on the location of the molecule being probed on the AI2O3 
substrate [17] that has a complex microscopic structure. 26|] In order to account for all of 
the data, we find that it is necessary to include the possibility of breaking of the fourfold 
symmetry of the Hamiltonian of the isolated molecule in the model. This is done phe- 
nomenologically in two different ways: In one of these (Approach A) it is assumed that 
the symmetry is broken by the interaction of the molecule with the complex substrate. In 
the other (Approach B) it is assumed that the symmetry is broken through the application 
of bias between the STM tip and substrate. We find that both approaches are generally 
successful but that Approach A is able to better account for one of the features of the ex- 
perimental data than Approach B. At present, since only one experiment of this kind is 
available, it is difficult to judge whether this difference between the two approaches confers 
a substantial advantage to one of them over the other. The approaches do, however, offer 
different predictions for experiments that have not yet been carried out. 

We find that photon emission is sensitive to the details of the molecule-substrate coupling, 
consistent with experimental data and the local-electrode interpretation of the experimental 
system. We also present calculated I-V characteristics for various coupling configurations 
and examine the relationship between the features found in the I-V characteristics and the 
occurrence and nature of the luminescence emitted by the molecule. Experimentally, photon 
emission was found to occur when there are two peaks in dl/dV. We find that for some 
coupling configurations, photon emission is predicted and the characteristic two-peak curve 
is obtained. For another configuration, only one peak in dl/dV is obtained and photon yield 
is very low. This is also in good qualitative agreement with experiment. Finally, we present a 
case of very weak molecule-substrate coupling, that has not yet been achieved experimentally, 
in which relatively high quantum efficiencies are predicted for photon emission. 

The organization of this article is as follows: In Section [III we describe our model and our 
method of solution. In Section UTTl we present our results, compare them with the experimen- 
tal data of Qiu, Nazin and Ho [17] and offer some predictions that may be tested in future 
experiments. In Section [IV] we present a concise summary of the aspects of the experimental 
data that our theory has been able to explain and of the physical mechanisms that we have 
identified as being responsible for them. We also comment further on the significance of the 
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present work for the fields of single-molecule electronics and optoelectronics. 



II. THE MODEL 



22 



The present model is a generalization of the simpler models presented in Refs. Il9| and 



. In the present model, the tip and substrate are represented by a tip electrode (probe) 
ectrodes (contacts), each modelled as one-dimensional tight-binding chains. 



and substrate e^ 
Unlike in Refs. 



19 



and 



22l where the formalism only allows single substrate contacts, in the 
formalism presented here an arbitrary number of substrate contacts are allowed. (We con- 
sider cases of 4 substrate contacts in the Results section of this article.) The roughly planar 
molecule lies on the substrate and is positioned between the tip and substrate electrodes, so 
that it mediates electron flow between the tip and substrate. The electronic model Hamilto- 
nian for this system can be divided into three parts, H = Heiectrodes + Hmoiecuie + W , where 
W is the interaction Hamiltonian between the electrodes and the molecule. Generalizing 



the Hamiltonian of Ref. 
electrodes is given by 



22 



to allow multiple substrate contacts, the Hamiltonian for the 



-1 



Heiectrodes= ^ e\n) {n\ + (3 {\n) {u - 1\ + \n ~ 1) {n\) 



+ ^ ^ + P{\n, i) {n + 1, i\ + \n + 1, i) {n, i\) , (1) 

i=l n=l 

where e are the site energies for the electrodes, (3 is the hopping amplitude between nearest- 
neighbor electrode atoms 28|, and \n) and \n,i) represent orbitals at site n of the tip probe 
and site n of the ith substrate contact, respectively. We take the electrochemical potentials 
of the tip and substrate electrodes to be /i^ = Ep + eVuas/'^ and iis = Ep — eVhiasf^, where 
Vhias is the bias voltage applied between them and Ep is their common Fermi level at zero 
bias. The Hamiltonian of the molecule may be expressed as 



molecule 



22 



(2) 

molecular 



where ej is the energy of the jth molecular orbital {\(f)j))- Unlike in Ref. 
orbital energies are allowed to shift in response to an applied bias voltage. Our treatment 
of the effect of bias voltage on orbital energies is described in Sec. IID. The interaction 
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Hamiltonian between the electrodes and molecule is given by 

m 

j 1=1 
where W-ij, Wj-i, and W(^i^i)j are the hopping amplitude matrix elements between 

the electrodes and the various molecular orbitals \4>j)- 

Electrons propagate in the form of Bloch waves through each electrode toward the 
molecule, and may undergo transmission or reflection when they encounter the molecule, 
contributing to the occupation of molecular orbitals in the process. Wave functions of elec- 
trons incoming from the tip probe are of the form 

— 1 m oo 

71=— oo i=l n=l j 

where d is the lattice spacing, ti are the transmission coefficients into the different substrate 
contacts, and r is the reflection coefficient. 



A. Solving the system 

In order to calculate molecular-based photon emission and I-V characteristics, it is neces- 
sary to evaluate the molecular orbital coefficients and transmission amplitudes for incoming 
electrons. This may be done by solving a Lippmann-Schwinger equation for this system, in 

but generalized to multiple substrate contacts: 



a similar fashion to Ref. 
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\(f>o) + Go{E)W\ilj), (5) 



where Go{E) = {E — {Heiectrodes + Hmoiecuie) + is the Green function for the decoupled 
system (without W), and |0o) is the eigenstate of an electron in the decoupled tip probe 
(or, more generally, the incoming electrode). Gq{E) may be separated into the decoupled 
components: the tip and substrate electrodes, and the molecule. For each electrode, 

^electrode ^ IM^)) {Mk)\ /„s 

° ^ E-{e + 2(3cos{kd)) ^ ^ 



k 



where d is the lattice spacing and e + 2(3cos{kd) is the energy of an electrode electron with 
wave vector k. ([^eZerfrode j^g^y g^jgQ expressed in an atomic orbital basis: 



oo oo 



^dectrode ^ ^ ^ (G'g'-*™^^)„ |n) (m | , (7) 



n=l m=l 



whose matrix elements {G\ 



electrode\ 
J 



are known analytically. 29|] For the molecule 



\4>j){4>j\ 

E - 



(8) 



For an electron incoming 
for the coefficients of \'ih): 



Tom the tip probe, this yields the following set of linear equations 



29| 



^-1 = {<t>o)-i + 



(Gf),(iy,,_iV^_i + W^,,(i,,,)V^i,) 



(9) 
(10) 

(11) 



where ip^i = ipi^i = and (0o)-i = (~l|0o)- The transmission probability for 

an electron incoming from the tip probe is given by T = J2T=i vl^iP-l28l Tip probe electrons 
between fix and in energy contribute to the electric current through the molecule. Using 
the Landauer theory, 30|] an expression for the current is obtained: 



2e /"'^^ 
/ = -/ T{E,VMas)dE. 



(12) 



The dependence of T on Vbias is due to shifting molecular orbital energies (described in Sec. 
IID). 



B. Photon emission 

Photon emission from the molecule can be understood in terms of allowed electronic 
transitions from a molecular orbital to one with a lower energy. To calculate emission 



spectra, as in Ref. 



19 



we use the expression for the spontaneous emission rate of a system 
emitting photons into empty space, using Fermi's Golden Rule. 3l| The emission rate is 
given by 



2, ,3 



3hc^ 



■mm)? 



(13) 



where ipi and ijjf represent initial and final states, and hw is their difference in energy. 



Unlike in Ref. 



19 



where photon emission is calculated for the idealized case of a two- 



orbital molecule, here we calculate photon emission for a system involving a real molecule 
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with multiple molecular orbitals. In order to do this, we consider emission only from the 
molecule itself. The rate is therefore approximated by 

R{ki,uj) = -^|^|c,v,^|2|c^.^|2|(0^,|x|0.)|2, (14) 
jd' 

where i and / label initial and final states. The relevant transition dipole moments (0j'|x|^j) 
are calculated by performing an extended Hiickel dipole analysis of the molecular orbitals. 32| 



To calculate the emission rate as a function of photon energy, we generalize the procedure 



presented in Ref. |19|. We must consider all electron states of the system, incoming from both 
the tip probe and each of the substrate contacts. Each electron state consists of an incoming 
wave, transmitted wave, reflected wave, and an amplitude on the molecule. See Fig. [T]for a 
schematic illustration. We assume here (and throughout the article) the positive bias case, 
with fix > fJ'S- Since we assume the temperature to be K, all states with incoming waves 
from a given electrode are occupied up to the electrochemical potential of that electrode. 
For a transition to occur, tpf must be an unoccupied state, and it must be lower in energy 
than ipi. Therefore, we consider transitions from occupied initial states (below /it) that are 
incoming from the tip probe, to unoccupied final states within the electrochemical potential 
window (above ^s) that are incoming from one of the substrate contacts. After normalizing 
the wave functions and converting the sum over k states (and spin) into an integral over 
energy, an expression for the photon emission spectrum (for a given bias voltage) is obtained: 

where Ei and ki are the initial energy and wave vector of an electron incoming from the tip 
probe, and lj is the frequency of the photon emitted. 



C. Electronic structure of Zn-etioporphyrin at zero bias 

The electronic structure of the Zn-etioporphyrin molecule was computed using the ex- 



tended Hiickel model. [25|| Within this model, the energy of the highest occupied molecular 
orbital (HOMO) is -11.5 eV, and the energy of the LUMO is -10.0 eV. For a charge neutral 
molecule at equilibrium, weakly coupled to the electrodes, the Fermi level of the electrodes 
at zero bias is expected to be located between the molecular HOMO and LUMO levels. 
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However, the precise location of the Fermi level is a difficult problem in molecular elec- 
tronics, with differing theoretical approaches yielding differing results. In STM experiments 
on Zn-etioporphyrin, the appearance of a low-bias dl/dV peak for some positions of the 
molecule above the substrate implies a Fermi level that is close to either the HOMO or 
LUMO. Since these are STM experiments, it is likely that the low-bias peak corresponds to 
an orbital entering the Fermi energy window by crossing ht (the electrochemical potential 
of the STM tip) rather than fis, due to the weaker coupling of the orbital to the tip than to 
the substrate. Since these experiments were performed at positive bias (electron flow from 
tip to substrate), it is therefore likely that the Fermi level is close to the LUMO and not the 
HOMO. 

The precise location of the Fermi level of the electrodes below the LUMO energy is likely 
to depend on the local geometry of the Zn-etioporphyrin/Al2O3/NiAl(110) interface: A 
work function study 33| of AI2O3 on NiAl(llO) has found that the formation of an ultra- 
thin AI2O3 layer on NiAl(llO) decreases the work function of the substrate by about 0.8 
eV, with a strong dependence on the oxide layer structure and thickness. It is therefore 
reasonable to assume that, for the experiments by Qiu, et al., ITj] different locations on the 
Al2O3/NiAl(110) substrate have different local work functions, with differences on the or- 
der of a few tenths of an eV. Due to these differences, variations in the common zero-bias 
Fermi energy of the tip and substrate (relative to the vacuum and also to the energies of 
the molecular orbitals) are likely to occur. Our calculations show that in most cases the 
overall qualitative picture is not sensitive to the precise location of the Fermi level below the 
LUMO energy. Therefore, for a qualitatively reasonable analysis of this system, we choose 
a zero-bias Fermi level of -10.1 eV. We justify this reasoning more explicitly in Sec. IIIB(3), 
where we compare our results for this Fermi level with results obtained assuming a zero-bias 
Fermi level of -10.3 eV. 



D. Molecular orbital energy-level dependences 

In order to realistically model photon emission and electric current as a function of bias 
voltage, it is necessary to consider the effects of bias voltage on molecular orbital energies. 
When a bias voltage is applied, an electric field is created between the tip and substrate, 
which may result in some charging of the molecule. If this occurs, the charging causes 
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an electrostatic shift of the molecular energy levels that in turn severely limits the actual 



charging that takes place. 



Ref. 



341 ] Generalizing the minimal charging approximation presented in 



19|, we phenomenologically approximate the shift of the molecular levels in response to 



the applied bias by adjusting ej for each molecular orbital so as to maintain the net charge 
that the molecule has at zero bias. The net electronic charge is calculated by summing 
over all occupied electron states (including spin) incoming from each electrode. This sum is 
converted into an integral, and an expression for the charge is obtained: 

^-2^7 2:.^^— (^^^+ y l^^p—(kdydE), (16) 

j ^ \ J contacts j ^' V / 

where the molecular orbital energies ej (and therefore cj) change with Vbias in such a way 
that Q = constant. 



1. Approach A 

One approach to treating the bias dependence of molecular orbital energies (which we 
will call Approach A) is to assume an equal bias dependence for the shifts in energy of each 
molecular orbital. This simple approach to charging yields physically reasonable behaviour 
of the molecular orbital energy levels with bias. However, by itself, it is insufficient to 
explain many of the experimentally observed STM I-V characteristics and photon emission 
results for the molecule Zn-etioporphyrin on Al2O3/NiAl(110). This may be because Zn- 
etioporphyrin has a two-fold degenerate lowest unoccupied molecular orbital (LUMO) that 
is likely to lose its degeneracy when the molecule is placed on a region of the complex 
surface where the molecule-substrate interaction is not fourfold-symmetric. After including 
such a substrate-dependent splitting in the zero-bias electronic structure of the LUMO, this 
approach yields interesting results that are consistent with the experimental data. 



2. Approach B 

Another approach (Approach B) is to consider the bias dependence of the different orbital 
energies in a slightly more complex way. Since Zn-etioporphyrin is a planar molecule and all 
of the relevant orbitals except for the LUMO have fourfold symmetry, the electric field from 
the STM tip affects each of the fourfold symmetric orbitals similarly, and we adjust their 
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energies by equal amounts a. The LUMO, however, consists of two degenerate orbitals with 
twofold symmetry. Depending on the position of the STM tip above the molecule, as bias 
voltage is applied this may result in a stronger electric field effect on the energy of one of 
the LUMO orbitals, and a weaker effect on the other orbital. Therefore, for cases where the 
tip probe is positioned above a region of the molecule with a high amplitude for one LUMO 
orbital and a low amplitude for the other, instead of adjusting the LUMO energies by equal 
amounts a we adjust the LUMO energies by amounts of 71 (> a) and 72 (< a) respectively. 
Within the present model, the quantities a, 71 and 72 all change with Vuas such that the 
total molecular charge Q remains constant. These quantities depend on the electrostatic 
geometry of the system. Therefore, for all values of Vuas-, the ratios a : 71 : 72 are kept the 
same, consistent with the linearity of electrostatics. With this phenomenological approach 
to charging, unlike in Approach A, we do not assume there to be any zero-bias splitting of 
the LUMO degeneracy. 

In the remainder of this paper we will present photon emission results and current- voltage 
(I-V) characteristics for Zn-etioporphyrin, calculated based on the above model for both 
Approach A and Approach B. We will show how photon emission is sensitive to details of 
the molecule-substrate coupling, and explore the relationship between photon emission and 
I-V curve features. In addition, we will demonstrate how our model can account for many 
previously unexplained features of the experimental data for this system. 



III. RESULTS 



We present results for Zn(II)-etioporphyrin I, coupled to a tip probe and 4 substrate 
contacts that we represent for simplicity by Cu s orbitals. The geometrical structure of the 
molecule has been calculated using density-functional theory. 35|] The molecule is mainly 
planar and oriented approximately parallel to the substrate, but contains four out-of-plane 
ethyl groups. 



A. Strong fourfold-symmetric molecule-substrate coupling 

We first consider a case where there is strong electronic molecule-substrate coupling rela- 
tive to the coupling between the molecule and the STM tip, and where the molecule-substrate 
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interaction is fourfold-symmetric. By 'strong coupling' we mean that the Hamiltonian ma- 
trix elements Weiectrodej between the relevant molecular orbitals and substrate contacts are 
about an order of magnitude greater than between the molecular orbitals and tip probe. 
It has been previously shown that the out-of-plane ethyl groups of the molecule are likely 
locations of dominant molecule-substrate coupling. 22| Therefore, four local substrate con- 
tacts (S1-S4) are positioned below the ethyl groups of the molecule, as shown in Fig. [2l|36| 
For Approach A (described in Sec. IID(l)), in this case we assume there is no splitting 
of the LUMO degeneracy, consistent with the fourfold symmetry of the molecule-substrate 
coupling. The tip probe is positioned (see Fig. [2]) above the molecule in a lateral region that 
has been shown to be part of the observed high-transmission lobe pattern for the STM tip 



above Zn-etioporphyrin. 



17 



22l | For this position of the tip probe (and any position corre- 



sponding to an experimentally observed high-transmission lobe) the tip probe has a stronger 
electrostatic coupling to one of the degenerate twofold symmetric LUMOs than to the other, 
and an intermediate coupling to all other relevant orbitals. [The difference between electro- 
static and electronic coupling should be noted: Electrostatic coupling refers to the change 
in the electrostatic potential that an electron in a molecular orbital feels due to the applied 
bias voltage, whereas electronic coupling refers to the Hamiltonian matrix element Weiectrodej 
between an electrode and a molecular orbital. In the rest of this article, these terms will be 
frequently used.] Therefore, for Approach B (discussed in Sec. I1D(2)), in order to model 
the shift of molecular orbital energies due to electrostatic effects in a phenomenological, 
quahtatively reasonable way, we assume the ratio a : 71 : 72 (discussed in Sec. IID) to be 
3:4:2. Here, 71 corresponds to the LUMO orbital that has stronger electrostatic coupling, 
and 72 to the orbital that has weaker electrostatic coupling to the tip. Results presented 
throughout this article are not sensitive to the precise values chosen for this ratio. 



1. Approach A 

For this strong substrate coupling case, with Approach A, photon emission is computed 
to be very weak. (In Sec. 1IIB,C, cases will be presented where the photon yield is more 
than an order of magnitude greater.) This weak emission result is consistent with the 
quenching of emission due to asymmetric coupling of the molecule to the tip and substrate 
observed experimentally jisl. and predicted for the general case of current-carrying molec- 
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ular wires [1 9(1 ■ The quenching of photon emission due to asymmetry of the electronic couphng 
can be understood physically as follows: Looking at FiglH in a highly asymmetric system 
where the tip-molecule coupling is much weaker than the molecule-substrate coupling, elec- 
trons incoming from the tip have relatively low amplitudes for entering the molecule, and 
high amplitudes for exiting into the substrate. There is therefore a low amplitude Cj^i for 
an electron in its initial state to be on a molecular orbital (even if the orbital is inside the 
Fermi energy window and close in energy to the energy of the electron), resulting in a low 
photon emission rate (see Eq jMj) . 

A further possible consideration is the molecular orbital amplitude cjj of an electron in 
its final state. If no allowed molecular orbitals are available to receive transitions (ie. inside 
the Fermi energy window of the system), cjj will be small for all possible final states and 
emission will be further quenched. As we will now show, for the strong fourfold-symmetric 
coupling situation we consider here, this in fact is the case. 

This further quenching, as well as the calculated current- voltage (I-V) curve for this case 
shown in Fig. [3^, can be understood by studying how the molecular orbital energies shift 
with bias voltage (see Fig. [3b): The LUMO (assumed to be degenerate in this case), becomes 
partially (slightly) occupied at low bias as the tip electrochemical potential (fir) approaches 
its energy. This causes an electrostatic shift of the molecular orbital energies (discussed in 
Sec. IID). The LUMO then shifts upwards in energy, following ht, so that the net charge 
on the molecule is maintained. 

The result is the approximately linear I-V curve at low bias in Fig. [3|l, with electron 
flow being mediated by the tails of the HOMO and the LUMO. There is a slight low-bias 
dl/dV feature due to /i^ approaching the LUMO energy. At about 1.3 V, the slope of 
the I-V curve begins to increase, resulting in a peak in dl/dV. The reason for this is as 
follows: The HOMO begins to become partially (slightly) unoccupied, even though it is still 
below the substrate electrochemical potential (/is). This is because the molecule-substrate 
contact couplings are strong compared to the molecule-tip coupling, so the substrate has a 
much stronger effect on the orbital occupations than the tip, and the high-energy tail of the 
HOMO begins to depopulate. HOMO electrons inside the Fermi energy window contribute 
to current flow into the substrate, increasing the slope of the I-V curve. The orbital energies 
are affected slightly, with the LUMO shifting slightly lower relative to fix (but not visibly 
in Fig. [3b), such that the net charge on the molecule is maintained. The slight downward 
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shift of the LUMO energy further increases the slope of the I-V curve. Here, electric current 
is very sensitive to such a shift, due to the LUMO's energy being very close to fiT- At 
about 1.4 V, the LUMO fully enters the Fermi energy window, in the process becoming only 
slightly occupied due to the much weaker coupling of the molecule to the tip (electron source 
electrode) than to the substrate (drain). At this point both the HOMO and LUMO orbital 
energies shift downwards, in such a way that the charge on the molecule remains constant (ie. 
the HOMO energy follows fis- The HOMO energy remains below /i^, resulting in quenched 
photon emission. The I-V curve flattens, since no orbitals are entering or approaching the 

energy window between tip and substrate Fermi energies. ^ ^ 

We now compare this result with experimental results obtained by Qiu et al.[17] for the 
STM/Zn-etioporphyrin/Al2O3/NiAl(110) system. In these experiments, depending on the 
location of the molecule on the substrate, the molecule either luminesced or did not, with 
different dl/dV curves obtained for luminescent and non-luminescent cases. See Fig.HJfor the 
reproduced experimental curves. Here, curve A and curve B are representative of molecules 
that were found to luminesce. Molecules with current-voltage curves C-F did not exhibit 
observable luminescence. Experimentally, molecules that did not luminesce were found to 
have only one dl/dV peak, usually at around 1.4 V. This is in good qualitative agreement 
with the model result presented here, using Approach A, which shows only one significant 
dl/dV peak that occurs at 1.4 V in Fig. [3^, and very weak photon emission, that is likely 
not experimentally detectable. 

2. Approach B 

For the case of strong fourfold-symmetric molecule-substrate coupling. Approach B (dis- 
cussed in Sec.IID) yields I-V results shown in Fig. |3b that are qualitatively similar to those 
in Fig. [3t that were obtained using Approach A. Photon emission is also computed to be 
very weak, for the same reasons as with Approach A. 

With Approach B, the LUMO with the weaker electrostatic coupling to the tip (which 
we will refer to as LUM02) enters the Fermi energy window at low bias (see Fig. [3ll), but 
contributes very little to the electric current (see Fig. [3b), due to the very weak LUM02-tip 
probe electronic coupling. The LUM02 remains almost completely unoccupied because of 
the asymmetry of the LUM02-tip and LUM02-substrate couplings. As approaches the 
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energy of the more strongly electrostatically and electronically coupled LUMO (LUMOl), 
however, the LUMOl becomes partially (slightly) occupied and shifts in energy, following /i^, 
so that the net charge on the molecule is maintained. The result is again an approximately 
linear I-V curve, with electron flow being mediated by the tails of the HOMO and the 
LUMOl. 

At about 1.5 V, the HOMO begins to become partially (slightly) unoccupied, similarly to 
Approach A, increasing the slope of the I-V curve. The LUMOl shifts slightly lower relative 
to /iT, such that the net charge on the molecule is maintained. This further increases the 
slope of the I-V curve. At 1.6 V, the LUMOl fully enters the Fermi energy window, in the 
process becoming only slightly occupied due to the asymmetry of the coupling. As with 
Approach A, the orbital energies then shift downwards, in such a way that the charge on 
the molecule remains constant. The HOMO energy remains below /i^, resulting in quenched 
photon emission, and the I-V curve flattens. 

For this case of strong molecule-substrate coupling using Approach B, there is found to be 
only one significant dl/dV peak (at 1.6 V) and very weak photon emission. As for Approach 
A, this compares well with the experimental non- luminescent cases (see Fig. HP-F) , where 
one dl/dV peak is observed (at about 1.4 V). 



B. Localized strong coupling 

Next, we consider the case where there is strong electronic coupling between the molecule 



and only one of the four substrate contacts. It has been suggested 



22| | that this type of 



electrode configuration is a likely possibility for the common experimental case of Fig.2B 
in the article by Qiu et al.[l3] Significant molecular electroluminesce was observed for this 
experimental case. 

The electrode configuration that we consider is similar to Sec.IIIA (see Fig. [2]); however, 
in this case the substrate contacts 5*1, 5*2 and 5*3 are moderately coupled to the molecule 
(couphng less than an order of magnitude greater than the coupling to the tip probe), and 6*4 



is strongly coupled. 38|] The tip probe is positioned in the same lateral region as for Sec.IIIA, 
and again with greater electrostatic coupling to one LUMO (LUMOl) relative to the other 
LUMO (LUM02). It should also be noted that, due to the twofold symmetry of the LUMO, 
the strongly coupled substrate contact is electronically strongly coupled to only one of the 
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LUMOs (LUM02, in this case) and not the other (LUMOl). 



1. Approach A 

With Approach A, since the molecule-substrate interaction is in this case not fourfold- 
symmetric, there is a splitting in the zero-bias degeneracy of the LUMO.[39] 

For this case, significant photon emission is computed to occur. Fig. shows the cal- 
culated emission spectrum at high bias {Vbias = 1-94 V). The spectrum corresponds to 
HOMO-LUMOl (1.94 eV peak) and H0M0-LUM02 (1.44 eV peak) transitions. The cal- 
culated I-V curve for this case, shown in Fig. [St, has a low-bias dl/dV peak and a high-bias 
dl/dV peak. 

To understand the calculated photon emission spectra and I-V curves for this case, it 
is necessary to pay close attention to the details of the coupling of the various molecular 
orbitals to the electrodes. Looking at FigJS^, at low bias the LUM02 enters the Fermi energy 
window, remaining almost completely unoccupied due to the strongly asymmetric coupling 
of the LUM02 to the tip and substrate. In this case, however, the LUM02 contribution to 
the electric current is not negligible. Current flow mediated by the LUM02 is not drowned 
out by current flow mediated through the tails of the LUMOl or the HOMO, since in this 
case the electronic coupling of the substrate is strongest to the LUM02. This creates the 
low-bias dl/dV peak seen in Fig. [5b. 

The LUM02 energy follows fix up to 0.2 V (see Fig. [5b). In this case, the substrate 
contacts have a large influence on the occupation of the LUM02 even though the LUM02 
is well above fis, because the coupling between the substrate and LUM02 is much stronger 
than between the tip and LUM02. Thus, from 0.2 V to 0.6 V the LUM02 tracks fxs and the 
I-V curve (Fig. \&p) is fiat. At 0.6 V, the LUMOl approaches /i^ and begins to populate. In 
response, the energies of the orbitals rise such that no charging takes place. The tip probe 
has a large influence on the occupation of the LUMOl, because the coupling between the 
tip/substrate and LUMOl is not highly asymmetric. From 0.6 V to 1.9 V, The LUMOl 
and the tail of the HOMO are the dominant sources of rising current. 

The HOMO reaches fis at Vbias = 1.9 V, causing an electrostatic shift in energy of the 
orbitals downwards, so that the LUMOl enters the Fermi energy window and populates sig- 
nificantly. The HOMO reaches fis and depopulates by an equal amount. There is a resulting 
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sharp increase in current, as both the HOMO and LUMOl mediate electron transmission 
from tip to substrate. 

Close inspection of Fig. and Fig. [3t shows that, in this case, the HOMO energy comes 
up to /i5 (within the resolution of the figure) whereas for the case of Sec. IIIA (Approach 
A), the HOMO energy only approaches fis- Hence, the depopulation of the HOMO is much 
greater in this case than in Sec. IIIA. Greater depopulation of the HOMO occurs, along with 
greater population of the LUMOl as it enters the Fermi energy window, due to the lack 
of strong coupling asymmetry between the LUMOl and the tip/substrate electrodes in this 
case. In this way, the zero-bias charge is maintained. Therefore, above 1.9 V, HOMO- 
mediated electronic states are available to receive transitions from LUMOl and LUM02- 
mediated states, resulting in photon emission. Since there is a stronger coupling asymmetry 
between the LUM02 and the tip/substrate electrodes than between the LUMOl and those 
electrodes, LUM02-H0M0 photon emission is weaker than LUMOl-HOMO emission (see 
Fig. [SK), as explained in Sec. IIIA. 

The onset of photon emission in this case occurs as the HOMO becomes partially unoc- 
cupied, at about 1.9 V. Fig. Eti shows the onset of photon emission, at the spectrum peak 
corresponding to H0M0-LUM02 transitions. Notice that the peak photon energy (^ 1.43- 
1.44 eV) is significantly less than the Fermi gap energy (^ 1.93-1.94 eV). This is because the 
LUM02 is deep inside the Fermi energy window at the onset voltage (see Fig. [5t). The cal- 
culated photon emission peak due to HOMO-LUMOl transitions has the same onset voltage 
as the H0M0-LUM02 emission peak. For this transition, however, photon energy is peaked 
close to the Fermi gap energy (1.9 eV) because the LUMOl and HOMO have energies close 
to /iT and /i5 respectively at the onset voltage. 

Comparing results for this luminescent case to experiment, the similarities are striking. 
Experimentally, molecules that luminesced commonly had a small dl/dV peak at 0.2 V and 
a larger peak at around 2.0 V (see Fig. |4]A,B). This is in excellent qualitative agreement 
with FiglSb, where we see a small dl/dV peak at 0.2 V and a larger peak at about 1.9 V. 
Furthermore, experimental results[l3] (reproduced here in Fig. [Tj) show the onset of photon 
emission occurring most commonly at about 2.2 V, but with a photon energy peak in the 
spectrum about 0.5 eV below the corresponding Fermi gap energy of 2.2 eV. This is in 
good agreement with Fig. [6l where at onset we find an emission peak (corresponding to 
the H0M0-LUM02 transition) significantly below the Fermi gap energy. Also, comparing 
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Fig. [6] with emission onset spectra for the most common experimental case (Fig. [7]) we see 
very similar behaviour of the emission spectra tails: The high-energy tail has a sharp cutoff, 
while the low-energy tail does not. As bias voltage increases, the high-energy cutoff shifts 
upwards in energy by a similar amount. In our model, we also see this behaviour, because 
the Fermi energy of the substrate provides a sharp energy cutoff below which there are no 
available final states for a transition. This cutoff reduces the extent of the high-energy tails. 
There is no such cutoff reducing the extent of the low-energy tails. 

Notice also that, experimentally, there is a shift in the position of the high-bias dl/dV 
peak, depending on whether photon emission is observed: In Fig. H] a peak is observed at 
1.4 V for non-luminescent cases (C-F), and around 2.0 V for luminescent cases (A,B). We 
see the same sort of bias peak shift theoretically with Approach A: 1.4 V for Sec.IIIA (weak 
emission case) and 1.9 V for Sec.IIIB (strong emission case). In this way. Fig. [3^ is similar 
to Fig. HP-F, while Fig. [Sb is similar to Fig. |1K,B. 

Our model further predicts a stronger HOMO-LUMOl emission peak (the 1.94 eV peak 
in Fig. [5^) with the same onset voltage as the experimentally observed H0M0-LUM02 
emission peak, but with a higher peak photon energy, close to the Fermi gap energy = eVuas- 
The experimental photon spectra in Ref. ,17 do not extend to the photon energy energy 
range in which this emission peak is predicted to occur (2.2 eV photon energy for the 
experimental onset voltage of 2.2 V). An experimental study testing this prediction would 
be very desirable. 

2. Approach B 

With Approach B, as with Approach A, significant photon emission is computed to occur 
in this case. Fig. [5b shows the emission spectrum at high bias {Yuas = 1-94 V). The 
spectrum corresponds to HOMO-LUMOl (1.57 eV peak) and H0M0-LUM02 (1.30 eV 
peak) transitions. The I-V curve for this case, shown in Fig. [Sji, has a low-bias dl/dV peak 
and a high-bias dl/dV peak. 

Looking at FigJSf, at low bias the LUM02 enters the Fermi energy window. It remains 
almost completely unoccupied due to the strongly asymmetric coupling of the LUM02 to 
the tip and substrate, but as with Approach A it still contributes to the electric current. 
This results in the low-bias dl/dV peak seen in Fig. [5ll. At 0.2 V the energy of the LUMOl 
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reaches /it- This causes an electrostatic shift in the energy levels upwards, as shown in 
Fig. [5F. From 0.2 V to 1.6 V, the LUMOl and the tail of the HOMO are the dominant 
sources of rising current. The HOMO reaches fis at 1.6 V, causing an electrostatic shift in 
energy of the orbitals downwards, so that the LUMOl enters the Fermi energy window and 
populates significantly. Similarly to Approach A, the HOMO reaches and depopulates 
by an equal amount, resulting in a sharp increase in current. 

For the same reasons as were explained for Approach A, for Approach B at 1.6 V HOMO- 
mediated electronic states are available to receive transitions from LUMOl and LUM02- 
mediated states, resulting in photon emission. As with Approach A, LUM02-H0M0 photon 
emission is weaker than LUMOl-HOMO emission (see Fig. \5jp). Fig. EJ) shows the onset 
of photon emission, around Vbias = 1.6 V, at the spectrum peak corresponding to HOMO- 
LUM02 transitions. As with Approach A, the photon peak energy is significantly less than 
the Fermi gap energy. 

Qualitatively, I-V and photon emission results for Approach B are similar to results for 
Approach A, and compare similarly well to experiment. There is one exception: With Ap- 
proach B, there is no shift in the position of the high-bias dl/dV peak depending on whether 
or not photon emission is observed: A peak is predicted at 1.6 V for both luminescent and 
non-luminescent cases, due to the very similar molecular orbital energetics for luminescent 
(Fig. [5]F) and non-luminescent (Fig. [Sji) cases. Experimentally, there is a shift in the position 
of the dl/dV peak: around 1.4 V for the non-luminescent case and 2.0 V for the luminescent 
case (see Fig. H])). A similar shift is found theoretically with Approach A, due to the fact 
that the HOMO-LUMOl energy difference in the luminescent case (Fig. [5fe) is greater than 
the HOMO-LUMO energy difference in the non- luminescent case (Fig. 

The physical reason for this difference between Approach A and Approach B is that in 
Approach A the molecule-substrate coupling splits the LUMO degeneracy in the luminescent 
case but not in the non- luminescent case and this difference in electronic structure results 
in the different bias voltages at which the high bias peak in dl/dV occurs. By contrast, in 
Approach B the LUMO degeneracy is lifted in both the luminescent and non-luminescent 
cases, so that the electronic structure of the molecule and the bias voltage at which the high 
bias peak in dl/dV occurs is similar in the two cases. 
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3. Dependence on the zero-bias Fermi level 



Experimentally, different dl/dV curves are observed depending on the location of the 
molecule on the substrate (see Fig. Even among those molecules that luminesced 

(A and B), there are differences in dl/dV. It should be noted that, in our article, we have 
chosen a zero-bias Fermi level of -10.1 eV, and that variations in the Fermi level relative 
to the molecular levels at zero bias are likely, depending on the location of the molecule on 
the surface, due to local work function variations (discussed in Sec. IIC). The dashed line 
in Fig. [5li shows an I-V curve (using Approach B) for an alternate zero-bias value of Ep: 
-10.3 eV instead of -10.1 eV. Here, the low-bias dl/dV peak is at 0.5 eV, corresponding more 
closely to Fig. HJA. than Fig. |3J3. It is possible that the experimental differences in low-bias 
dl/dV peak locations in Qiu's Fig.2A and Fig.2B are due to different zero-bias Fermi levels 
caused by local work function variations on the surface. Other than the change in the low- 
bias dl/dV peak location, small changes in the Fermi level yield qualitatively similar I-V 
and photon emission results. Therefore, in the rest of this article, we have assumed a Fermi 
level of -10.1 eV. 



C. Weak fourfold-symmetric molecule-substrate coupling 

The final case we consider is weak molecule-substrate coupling, where the molecule- 
substrate interaction is fourfold-symmetric, along with stronger tip-molecule coupling than 
in the previous cases. In this case, the electronic molecule-substrate coupling is of the 



same order of magnitude as the tip-molecule coupling. [40|] This situation may be achieved 
experimentally by increasing the thickness of the oxide layer between the molecule and metal 
substrate by a modest amount, or by decreasing the tip-molecule distance. In our model, 
we both increase the molecule-substrate distance and decrease the tip-molecule distance. 



1. Approach A 

As in Sec.IIIA(l), we assume there is no splitting of the LUMO degeneracy. For this case, 
much more efficient photon emission is predicted to occur, with a photon yield two orders 
of magnitude higher than for Sec. IIIB. Fig. [8^ shows the emission spectrum at high bias 
{Vbias = 1-95 V). The peak in the spectrum corresponds to the HOMO-LUMO transition. 
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Fig. [Ht shows the I-V curve for this case. There is a high-bias dl/dV peak (at 1.45 V) and 
no low-bias peak. 

Looking at Fig. [S^, the molecular orbital energetics are similar to those for Sec.IIIA(l) 
(shown in Fig. [3b). Since no orbitals enter the Fermi energy window at low bias, there is 
no low-bias dl/dV peak. In this case, the I-V curve is quite flat up to about 1.4 V. At 
1.4 V, the LUMO fully enters the Fermi energy window, becoming partially occupied. The 
HOMO depopulates by an equal amount, and the orbitals electrostatically shift downwards 
in energy with ns- 

Since the tip has a much stronger effect on the LUMO occupation in this case than 
in Sec.IIIA,B, the degree of partial population of the LUMO, and partial depopulation of 
the HOMO, is much greater. This results in much greater quantum efficiency for photon 
emission. Unlike in Sec.IIIB, the initial onset voltage for photon emission due to HOMO- 
LUMO transitions matches the HOMO-LUMO emission peak energy. 

2. Approach B 

As with Approach A, with Approach B very strong photon emission is predicted to occur. 
Fig. [8]d shows the emission spectrum at high bias (Hias = 1-80 V). Unlike for Approach A, 
here there are two peaks in the spectrum, corresponding to HOMO-LUMO 1 and HOMO- 
LUM02 transitions. Fig. [8li shows the I-V curve for this case. There are two high-bias 
dl/dV peaks (at 1.2 V and 1.6 V) and no low-bias peak. 

These results can be understood by studying the behaviour of the molecular orbitals 
with applied bias voltage (see Fig. [8]F). At low bias, the Fermi energy window approaches 
the LUMOl and LUM02. Unlike the other cases (Sec.IIIA,B), in this case the LUM02 
coupling to tip and substrate is not strongly asymmetric, and electron states from the tip 
have a significant effect on the charge of the orbital; therefore, the LUM02 electrostatically 
shifts in energy with /i-p so that the zero-bias charge on the molecule is maintained. Since 
no orbitals enter the Fermi energy window at low bias, there is no low-bias dl/dV peak. At 
1.2 V, the HOMO energy reaches /x^, and the HOMO begins to depopulate. This causes 
an electrostatic shift in orbital energy downwards, and the LUM02 enters the Fermi energy 
window, creating a dl/dV peak at 1.2 V. At 1.6 V, the LUMOl enters the Fermi energy 
window, resulting in another dl/dV peak. (This increase in current is greater than the 
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increase at 1.2 V, because the LUMOl has stronger electronic couphng than the LUM02 to 
the tip probe.) The HOMO depopulates significantly further, with the LUMOl populating 
by an equal amount. (The resulting electrostatic deviation in orbital energies is too small to 
be visible in Fig. [H]F because the occupation of the HOMO is very sensitive to any deviation 
in energy away from /i^.) 

As with Approach A, the result is higher quantum efficiency for photon emission. Un- 
like in Sec.IIIB(2), the initial onset voltage for photon emission due to H0M0-LUM02 
transitions corresponds to the H0M0-LUM02 emission peak energy. The H0M0-LUM02 
emission peak increases further once the onset voltage corresponding to the HOMO-LUMOl 
emission peak is reached (due to the further depopulation of the HOMO). 

A signature of this relatively efficient photon emission regime, found with both Approach 
A and Approach B, is the lack of a low-bias dl/dV peak. This regime has yet to be realized in 
STM experiments; however, it is predicted that greatly enhanced quantum efficiency could 
be achieved by further weakening the coupling of the molecule to the metallic substrate, 
or by bringing the STM tip closer to the molecule. While in our model we both increase 
the molecule-substrate distance and decrease the tip-molecule distance, it may be more 
experimentally feasible to increase the thickness of the oxide layer without bringing the tip 
closer to the molecule. This would cause a reduced current through the molecule. For such 
an experimental situation, the relevant luminescence observation is not the absolute photon 
emission intensity, but the quantum efficiency, or photon yield (the number of photons given 
off per electron passing through the molecule). This is predicted to be greatly enhanced. 

D. Discussion of Results 

Both Approach A and Approach B yield results consistent with experiment. For the case 
where the molecule is strongly coupled to the substrate, very weak photon emission, along 
with only a single high-bias dl/dV peak, is found with both approaches. Experimentally, all 
molecules that did not luminesce had a single high-bias dl/dV peak signature and no low 
bias dl/dV peak. 

For the case where only a localized region of the molecule is strongly coupled to the 
substrate, both approaches yield much stronger photon emission than the first case. This is 
because, for a HOMO-LUMO transition, the relevant coupling asymmetry (between the tip- 



23 



LUMO and the HOMO-substrate) is greatly reduced. Two emission peaks were found, the 
lower-energy peak being significantly lower in energy at onset than the energy corresponding 
to the onset voltage. As well, in this case both a low-bias and high-bias dl/dV peak are 
found. This is consistent with experiment: In the experimental case where both low-bias 

urthermore. 
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and high-bias dl/dV peaks are observed, photon emission is also observed, 
there is additional evidence based on modelling of the molecular STM images 
experimental case corresponds to a localized region of strong coupling of the molecule to the 
substrate. 

One qualitative feature observed experimentally and found theoretically with Approach 
A is not found with Approach B: Experimentally, there is a shift in the position of the 
high-bias dl/dV peak, depending on whether or not photon emission is observed. This shift 
is predicted with Approach A but not with Approach B. 

There is additional experimental evidence in support of Approach A in the form of an 
observed zero-bias splitting in the LUMO degeneracy of a similar molecule (magnesium 
porphine) above the same Al2O3/NiAl(110) substrate. j4l| It should be noted that for this 
experiment, only MgP molecules with two- lobe STM images were chosen for detailed study, 
so the substrate-dependence of the zero-bias splitting is unknown. 

There is, however, a possible physical justification for Approach B. When a bias voltage is 
applied, the STM tip will electrostatically affect different molecular orbitals differently. The 
extent of these different effects is unknown. A simple electrostatic calculation, treating the 
tip/substrate as a point charge and a mirror image charge, suggests small differences (typi- 
cally on the order of lOOths of an eV) in the average potential for the LUMOl and LUM02 
orbitals. Thus, while the assumptions for Approach B may indeed be qualitatively correct, 
the degree to which the orbital energies of the LUMOl and LUM02 behave differently with 
bias is unknown and may be small. 

For the case of very weak molecule-substrate coupling, much higher quantum efficiencies 
for photon emission are predicted to occur. This regime has not yet been realized experi- 
mentally, and would be an intriguing avenue for further research. 
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IV. CONCLUSIONS 



The local-electrode framewor 
of experimental observations 



c presented in this article coherently explains a multitude 
17l | not previously theoretically studied, for the STM/Zn- 
etioporphyrin/Al2O3/NiAl(110) system. The following is a summary of these observations, 
with explanations based on our model results: 

(i) The observed molecular-based photon emission is due to transitions between the molec- 
ular LUMO, whose degeneracy has been split by molecule-substrate and/or molecule- 
STM tip interactions, and the molecular HOMO. 

(ii) For some cases, low-bias dl/dV peaks are observed experimentally (see Fig. |3]A.,B). 
Our model explains these as being due to a splitting of the LUMO degeneracy, with 
the lower-energy LUMO entering the Fermi energy window at low bias (see Sec. IIIB). 

(iii) For some cases, no low-bias dl/dV peak is experimentally observed (see Fig. HP-F). 
We find that this occurs because the molecule is too strongly coupled to the substrate, 
with the LUMO either not entering the Fermi energy window at low bias (Approach A, 
see Sec.IIIA(l)), or entering the window but contributing negligibly to the current due 
to very weak coupling of the molecule to the tip compared to the substrate (Approach 
B, see Sec. IIIA(2)). 

(iv) For cases with no low-bias peak, no photon emission is experimentally observed. This 
is due to strongly asymmetric tip/molecule and molecule/substrate couplings. In these 
cases, when a bias is applied, the HOMO stays almost fully occupied and the LUMO 
almost completely unoccupied (see Sec. IIIA(1,2)). 

(v) There is an experimentally observed difference in the position of the high-bias dl/dV 
peak, between cases where photon emission is and is not observed (see Fig. H]). This 
is explained with Approach A by a breaking of the LUMO degeneracy only in the 
luminescent case (see Sec. IIIA(l) and Sec. IIIB(l)). 

(vi) The experimental peak photon energy is about 0.5 eV below eVbias at emission onset 
(see Fig. [7j). This is due to splitting of the LUMO degeneracy, with the lower-energy 
LUMO being well inside the Fermi energy window as the energy of the HOMO ap- 
proaches the window. See Sec. IIIB(1,2). 
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(vii) The high-energy photon emission spectra tails are steeper than the low-energy tails 
(see Fig [7]). This is due to the substrate Fermi energy providing a sharp energy cutoff 
below which there are no available states to receive a transition (see Sec. IIIB(1,2)). 

(viii) There are significant differences in experimentally observed positions of dl/dV peaks 
(see Fig H]) depending on the position of the molecule on the substrate. These differ- 
ences are consistent with differing local zero-bias Fermi levels, due to local variations 
in the work function of the oxide-coated metal substrate (see Sec. IIIB(3)). 

Our model predicts an additional photon emission peak to be found, having a peak energy 
close to the bias voltage at emission onset, for the case of molecular-based photon emission 
presented in Sec. IIIB. Experiments testing this prediction would be of interest. 

We also predict that greatly enhanced quantum efficiency of photon emission could be 
achieved by further weakening the coupling of the molecule to the metallic substrate, or if 
possible by bringing the STM tip closer to the molecule (see Sec. IIIC). The spectrum for 
this greatly enhanced quantum efficiency could yield further clues to the relative merits of 
the two approaches studied in this article (see Sec.IID). For Approach A, one emission peak 
is predicted, and for Approach B, two peaks are predicted. 

Studying the STM/Zn-etioporphyrin/Al2O3/NiAl(110) system using the local-electrode 
theoretical framework presented in this paper has yielded a coherent explanation of a large 
body of experimental results for this system. Using this framework, we are able to gain a 
much greater understanding of single molecule electroluminescence. This is an important 
step towards the development of the emerging field of single molecule optoelectronics. We 
hope that this work inspires further experimental and theoretical research in this promising 
new field. 

While the present theory relies heavily on phenomenology, it has allowed us to construct 
energy level diagrams of the evolution of the molecular HOMO and LUMO orbitals and 
of the electrochemical potentials of the electrodes as a function of applied bias that are 
physically reasonable and are consistent with both the experimentally observed current- 
voltage characteristics and the experimental electroluminescence data. Thus the present 
work can also be viewed as a quantitative interpretation of the experimental data that is 
unique in that it satisfies more demanding experimental constraints than previous attempts 
to model experimental molecular electronic data that have focussed on experimental current- 
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voltage characteristics alone. Therefore, as well as contributing to a better understanding of 
single-molecule optoelectronics the present work provides much needed benchmarks for the 
development of accurate first principles theories of the evolution of the electronic structure 
of molecular nanowires under bias that do not yet exist at this time. 
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FIG. 1: (Color Online) A schematic energy level diagram of a transition from an occupied electron 
state (incoming from the tip probe on the left, shown in red) to an unoccupied electron state 
(incoming from a substrate contact on the right, shown in blue). A photon is created with energy 
huj, equal to the difference in energy between the two electron states. The dashed lines represent 
the molecular portions of the states. 

FIG. 2: (Color Online) The Zn(II)-etioporphyrin I molecule, showing substrate contacts Si, S2, S3, 
and S4 (open blue circles, into the page) and the tip probe (blue dot, out of page). Carbon atoms 
are red, nitrogen atoms are green, the zinc atom is yellow, and hydrogen atoms are white. 

FIG. 3: (Color Online) Strong fourfold-symmetric coupling between molecule and substrate: Elec- 
tric current, and molecular orbital energies as a function of bias voltage, (a) Approach A, I vs. 
Vbias- Red lines represent dl/dV. (b) Approach B, I vs. Vbias- (c) Approach A, molecular orbital 
energies (dashed lines represent tip and substrate electrochemical potentials), (d) Approach B, 
molecular orbital energies. 



FIG. 4: (Color Online) From Qiu et al.[l7t]. Reprinted with permission from AAAS. Experimental 
dl/dV curves for Zn-etioporphyrin/Al2O3/NiAl(110) obtained with the STM, for molecules at 
different locations on the substrate. (A-F) dl/dV curves representative of the various molecular 
images observed. The curve seen in B was most commonly observed (30% of the time). Molecular 
electroluminescence was observed for cases A and B but not for C-F. 

FIG. 5: (Color Online) Localized strong coupling: photon emission, electric current, and molecular 
orbital energies as a function of bias voltage, (a) Approach A, photon emission vs. Vbias- (b) 
Approach B, photon emission vs. Vbias- (c) Approach A, I vs. Vbias- Red lines represent dl/dV. 
(d) Approach B, I vs. Vbias- Dashed line represents Ep = —10.3 eV. (e) Approach A, molecular 
orbital energies, (f) Approach B, molecular orbital energies. 

FIG. 6: Localized strong coupling: Onset of photon emission, at the II0M0-LUM02 emission 
peak, (a) Approach A, emission rate vs. photon energy, for three different values of Vbias around 
the onset voltage, (b) Approach B, emission rate vs. photon energy. 
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FIG. 7: (Color Online) From Qiu et al.[l3]- Reprinted with permission from AAAS. Experimental 
photon emission spectra for molecules corresponding to FigjJjB, for various bias voltages around 
the onset voltage. 



FIG. 8: Weak fourfold-symmetric coupling: photon emission, electric current, and molecular orbital 
energies as a function of bias voltage, (a) Approach A, photon emission vs. Vbias- (b) Approach B, 
photon emission vs. Vhias- (c) Approach A, I vs. Vhias- (d) Approach B, I vs. Vhias- (e) Approach 
A, molecular orbital energies, (f) Approach B, molecular orbital energies. 
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